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COARSENING IN ONE DIMENSION: 
INVARIANT AND ASYMPTOTIC STATES 

EMANUEL A. LAZAR AND ROBIN PEMANTLE 


Abstract. We study a coarsening process of one-dimensional cell complexes. We show that if cell 
boundaries move with velocities proportional to the difference in size of neighboring cells, then the 
average cell size grows at a prescribed exponential rate and the Poisson distribution is precisely 
invariant for the distribution of the whole process, rescaled in space by its average growth rate. We 
present numerical evidence toward the following universality conjecture: starting from any finite 
mean stationary renewal process, the system when rescaled by e~ 2t converges to a Poisson point 
process. For a limited case, this makes precise what has been observed previously in experiments 
and simulations, and lays the foundation for a theory of universal asymptotic states of dynamical 
cell complexes. 
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1. Introduction 

Many physical systems can be abstracted as cell complexes whose geometry and topology change 
over time through deterministic evolution equations. Soap foams and polycrystallinc metals are two 
examples of such systems that evolve under a generalized curvature flow to reduce an energy asso¬ 
ciated with the co-dimension 1 cells. Over time, some cells grow while others shrink and disappear. 
This is sometimes called a coarsening process, because as cells disappear, the average size of the 
remaining cells monotonically increases [Lazllj . 

A remarkable self-similarity property has long been observed in both experiments and numerical 
simulations of these systems: if an arbitrary cell complex is allowed to evolve by fixed evolution equa¬ 
tions, then its scale-independent properties will converge, in a statistical sense, to a time-independent 
set of properties. For example, the distribution of normalized cell sizes converges to some fixed dis¬ 
tribution, even as the system itself continues to evolve |Lazlll iLMMSlTI [MLMS121 . This general 
behavior is sometimes referred to as statistical self-similarity |Mul861 IMV89] . A second remarkable 
universality property is also observed: these self-similar properties are largely independent of the 
initial conditions, and instead depend on the deterministic evolution dynamics. 

The primary aim of this paper is to give a rigorous analysis of a one-dimensional coarsening 
process in which a limit theorem can be proved. The second goal is to explore the degree to which 
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Figure 1. A one-dimensional cell complex with a few of the cell boundaries X{ labeled. 


the limiting behavior is independent of initial conditions. These will be carried out in the context 
of a particular one-dimensional system which can be viewed as one-dimensional curvature flow. 

We parametrize configurations by the locations x t of the boundaries between the 1-cells. The 
particular system we initially consider is one in which cell boundaries move with a velocity propor¬ 
tional to the difference of the sizes of two neighboring 1-cells. We let X = (x*) £ R be an ordered 
set of points on a one-dimensional manifold, as illustrated in Figure [T| The evolution of this system 
is governed by two rules: 

(1.1) ^ = -Xj_i + 2xi - x x+1 

(1.2) If Xi = Xj for i ^ j : we remove Xj and reindex the remaining points. 

Rule 2 prevents points from crossing one another. It is straightforward to see that 1-cells that are 
larger than their two neighbors grow, while those smaller than their neighbors shrink. 


Previous Work. Important work has been done in studying one-dimensional dynamical cell com¬ 
plexes, though these systems are often known by the names clustering or coagulation systems. Carr 
and Pego [CP92] considered a model in which at discrete time steps the shortest interval in a finite 
system is taken and joined with its neighbors. They use Laplace transforms to establish the existence 
of a self-similar solution. This system makes use of discrete time steps and “long-range” interaction. 
Derrida |Der951 IDer97| considered one-dimensional Ising and Potts models at zero temperature and 
considered the rates at which these systems coarsen. He showed that if the initial condition is ran¬ 
dom, the system enters a scaling regime where the average size of domains grows as t 1 / 2 . Nontrivial 

exponents that govern other properties of these coarsening systems are also calculated. Hunderi 

doc ■ 1 1 

and Ryum lHR96tIHRW79] considered a system in which —-L =-. They used 

dt Xi -|_i Xi Xi Xi —i 

numerical simulations to demonstrate that a size correlation develops between neighboring cells; this 
correlation calls into question the reliability of mean field theories I I1S92 . which tend to ignore 
size correlations between neighboring cells, in analyzing these systems. 

Mullins [Mul91| considered a slightly more general version of the system considered in this pa¬ 
per, and established some very preliminary results for these systems. Unfortunately, despite work 
of Mullins and of others, Mullins’ conclusion, “Very little is known, however, about the general 
conditions under which self-similarity is to be expected” remains true nearly twenty-five years later. 
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2. Results 

Our main result is the invariance of the Poisson distribution for cell boundaries under the evo¬ 
lution (11.11) (11.21) . We begin with an analysis of what happens with a finite number of such points, 
moving in a finite interval. For boundary conditions, here and throughout, we impose two extra 
particles that remain stationary at the endpoints 0 and 1. 

Theorem 2.1. Begin with n points chosen uniformly on [0,1]. The law after time t is that of Z 
points chosen uniformly and independently on [0,1] where Z is a binomial with parameters n and 


Because these points are moving around, colliding with each other and coalescing, we refer to 
them henceforth as particles. Some consequences follow immediately from Theorem 12. II properties 
of the binomial and Poisson distributions, and scaling invariance. Let V(X,L) denote the law of a 
Poisson process of intensity A times Lebesgue measure on [-L, L\; alternatively it can be described 
as the law of Z IID uniform points in [—L,L] when Z has Poisson distribution with mean A. 


Corollary 2.2. 


1. Starting from n particles uniformly distributed in [—L,L\, 
subset A C [—L,L\ after time t has binomial distribution 

„-2 1 l^-l 


Bin 


n, e 


2 L 


the number of particles in a given 


2. Starting from V{X,L), the law of the configuration after time t is a Poisson configuration with 
intensity Xe~ 2t times Lebesgue measure on [—L,L\. 

3. For M > L > 0, let C(M,L) be the law of the time t configuration restricted to {—L,L\, starting 
from V(X, M). Then the law C(M,L) does not depend on M. 


Next we turn to the infinite process. The first thing is to establish an existence result. This is not 
trivial because of the possibility of influence coming in from infinity in finite time and because at 
collision times, which are dense, the trajectory of every particle become non-analytic. In Section [7] 
we will define the space LI of infinite particle configurations and prove the following result. 

Theorem 2.3 (existence of a weak solution for the infinite system). Let LI be the space of continuous 
trajectories of infinite configurations defined in Section [?| For every A > 0 there exists a weak 
solution to (O) with initial condition Poisson of intensity X, in the sense of Definition \7.2[ Any 
such solution has the following properties. 

(i) The trajectory of each particle is always differentiable from the right. 

(ii) Excluding a measure zero set of initial configurations, the collision times will be distinct. 

(in) At all times that a particle is not involved in a collision, that particle’s position has a two-sided 
derivative. 
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The infinite model allows for a statement of Poisson invariance cleaner than the one in Corol¬ 
lary 12.21 

Theorem 2.4. Under any weak solution Pa, the time t law of the configuration is Poisson with 
intensity \e~ 2t . 


3. Evolution on a finite interval 

In order to define the process on a finite interval, we need to choose how the first and last particles 
evolve. As previously mentioned, rather than using periodic boundary conditions, we impose extra 
particles at the endpoints that never move. 

The space of configurations of n particles on [0,1] is the simplex 

A n := {x = (0, xi ,..., x„, 1) : 0 = xo < x\ < ■ ■ ■ < x n < x n +i = 1} ■ 

We include the degenerate simplex A 0 := {(0,1)}. The interior A° of A„ is the set of configurations 
with all particles distinct from each other and from the endpoints {0,1}. After a collision, two 
particles, say Xi and , get stuck together. The configuration is no longer in A°, rather it is 
now on the boundary dA n . The boundary is composed of faces of various dimensions defined below 
as images of open simplices under various embeddings. On each of these, the evolution satisfies a 
different rule. Rather than think of this as a single discontinuous rule on the closed simplex, we 
think of the particle as entering the (n — l)-simplex at the first collision time. 

To make this precise, we define a set of embeddings as follows. For integers n > 0 and k G [0, ra+l], 
define the map : A n —>• A n+ i by the formula 

! (0, 0, aii,..., x n , 1) k = 0 

(0, ny,..., x U f 1,1) k n T1 

(0, x\,.. .,x k ,x k , ...,x n ,l) 1 <k<n 

Define the projection map from A„ to the union of A£ for k < n by letting 7r(0, x\,... x n , 1) be the 
result of omitting repeated entries. Then tt o t n _ l fc (x) = x for each x e A"_j and each k G [0, n]. 
More generally, any composition of maps o o im,k m is inverted by n on Aj^. 

We now define an evolution on the configuration space A := (+J n A°, the disjoint union of n-point 
configuration spaces on [0,1]. The evolution is defined as a set of time t maps {$ t } on A and 
constructed by induction on the dimension n of the stratum in A. The base step of the induction is 
to define ^(x) e 0 on Aq. 
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Let A n denote the (n + 2) x (n + 2) matrix defined by 


0 0 0 0 ••• 0 

-1 2 -1 0 ••• 0 

0 -1 2 -1 ••• 0 


0 ••• 0 -1 2 -1 

0 ••• 0 0 0 0 


The differential equation x' = A n x defines a flow on A°. For x £ A°, let 


t(x) := inf {t : e tAn x £ 9A„} . 


We may now recursively define 

$ t (x) := 


e tAn x 


t < r(x) 
^t-r(x) (n(e rAn x)) t > r(x) 


Informally, flow by x' = A n x until you reach the boundary, then collapse whichever points need 
to be collapsed and continue inductively (run x' = A m x where m is the new number of particles, 
and so on). For all but a set of measure zero of initial points, each collision will reduce the number 
of particles by precisely 1; however, the flow is well defined even for simultaneous coalescences, in 
which case n drops dimension by more than 1. 


Remark. One can define the evolution directly on A n . At time t(x), let F be the unique face of 
A n to which x is interior. Let i be the composition of maps t mt k inverted by 7r on F\ the map l is 
unique even though it may be represented by a sequence of compositions in many ways. Instead of 
mapping x to the projection 7r(d> r ( x) (x)) of its time r image, continue the evolution via 

<f> T+s (x) := t ($ s (7r(e Tj4 "x)) . 

Informally, the evolution on A n keeps track not only of the coalesced system but of which original 
particles coalesced at what is now particle i. 


4. The time t pre-image of a finite configuration 

The forward flow goes until it hits the boundary, jumps down a dimension, and continues. The 
reverse flow goes inward away from the boundary, but unlike the forward flow, it is capable at any 
time of reinterpreting its position in A m as a position on dA m+ i , in which case it jumps up a 
dimension and begins to flow inward from the boundary of the bigger simplex. Not only that, but 
it can choose to jump up in any of m + 2 ways: any of the particles can split into two, or a particle 
can enter at 0 or 1. In this section we characterize all possible reverse trajectories. First we check 
that the flow, when not splitting, does indeed flow inward. 
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Lemma 4.1. For t > 0, the map x i —> e tj4 "x is one to one on and maps A n into its own 

interior. 

Proof: Letting yj = x 1+i — Xj, the ODE x' = — A„x induces an ODE y' = B n y where B n is the 
(n+ 1) x (n+ 1) matrix with ones on the first super- and sub-diagonal, and —2 on the main diagonal 
except that B\\ = B n+ i, n +\ = —1, not —2. Each ( B n y)j is a linear combination of coordinates of y 
with the only negative contribution coming from yj\ hence if y £ (R + )" and yj = 0 then {B n y)j > 0. 
By Nagumo’s Theorem Aub77] (or see lHar721 Theorem 1]), solutions to y' = B n y never leave the 
nonnegative orthant. It follows that the coordinates Xj never collide. □ 

Now fix integers n > to > 0 and a point y £ A(^, and define the following notation for the time 
t dimension n pre-image: 

$t -1 (y; n) := {x £ A n : y = $ t (x)} . 

We show that <E’ t ” 1 (y;n) has a one to one parametrization by the sequence of times at which the 
split occurs and the choice of which particle splits at each time. 

To define this formally, let T(f; k) be the fc-simplex of vectors t satisfying 0 < t\ < ■ ■ ■ < tk < t. 
Define 

S(n\m) := {(fci,..., k n - m ) : 0 < hi < n - i + 1} . 

As time passes backward through tk, one of the existing n — k particles splits; these choices are what 
is encoded by an element of S(n; m). 

Proposition 4.2. Fix integers m < n and a positive real t. Let k := n — m. There is a map p from 
A m x T(t;k) x S(n;k) to A n such that if x £ rj(y,t,s) then 

(i) $ t (x) = y 

(ii) For each 1 < j < k, the first entry time Tj of < I > .(x) into A„_j is equal to tj 
(Hi) For each 1 < j < fc the coalescence at time Tj occurs between coordinates Sj and Sj-\. If 
Tj = ■ ■ ■ = Tj +r , this is taken to mean that the r + 1 coalescences at this time are those 
identifying Sj+i with for i = 0,..., r. 

Proof: We construct the map rj explicitly. Run the clock backwards starting at time t, evolving 
backwards via e~ tAm . By Lemma l4Tl the point stays in A^. When the clock reaches tk, jump from 
A m to A m+ i via the embedding t m ,s k - Continue evolving backwards via e~ tAm+1 until the clock 
reaches time tk- i, then jump up to A m+ 2 via Continue in this manner, reaching A n at 

time t\ and evolving in A n by e~ tAn for a backward time t.\ until the clock says zero. 

With this construction of y, the consequences of the proposition are easily checked by induction 
on k. When k = 0 the vectors t and s are empty and the only fact to check is that e~ tA inverts e tA . 
For k > 1, evolving back in time by t — tk yields a point y' := e~ tAm y and a point y" := L m +i,s k y'■ 
Apply the induction hypothesis to x = ?y(y',t',s') where t' and s' drop that last coordinates of 
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t and s respectively to see that <f> tfc _(x) := lim., ttfc 3> s (x) = y". Thus, $t(x) = ($ ifc (x)) = 
$t-t k ( y') = y. One final coalescence occurs at time tk in coordinates Sk and Sk — 1. This completes 
the induction. □ 


5. Invariance theorem for finite configurations 

This section is devoted to the proof of Theorem 12.11 and its consequences. We reduce it to 
Lemma 15.21 below, prove the further Lemma 15.31 then prove Lemma 15.21 It will be helpful to define 
a certain continuous time Markov chain on the nonnegative integers. It is a pure death chain, that 
is, transitions from n are allowed only to n — 1. The rate at which n transitions to n — 1 is 2 n. If 
{X t : t > 0} is such a Markov chain, one interpretation is that X t is the number of particles alive 
at time t, where each particle independently dies at rate 2 (meaning, after an exponential random 
time with mean 1/2). Denote the transition probabilities for this chain by pt . Thus, pt(n,m) is the 
probability, starting with n particles, that precisely m are alive at time t. Saying it another way, 

(5.1) Pt(n, •) ~ Bin (n, e~ 2t ). 

Let g(t\, ..., tk) = g{n\ ti,... , tk) denote the density for the first k transitions starting from n. An 
explicit formula for g is 

(5.2) g(n-, h, ...,t k ) = 2 k n(n - 1) • • • (n - k + 1) exp (- ^ 2 (n-j + 1 )(tj - tj- 1 ) 

V i= i 

where to = 0 by convention. 

Normalized Lebesgue measure on A n is the probability measure p n whose density with respect 
to Lebesgue measure is n\. We first prove a special case of Theorem 12.II 

Lemma 5.1. Let B be any Borel subset of A n and let t be any positive real number. Then 

Hn{x e A n : r x > t and e ti4 x eB} = e~ 2nt p n (B). 

Consequently, after time t the probability of still having n points is e~ 2nt , and conditional on this 
the points are independently and uniformly distributed. 

PROOF: The divergence of the vector field F(x) = A n x on R n + 2 is equal to 2 n. It follows that the 
Jacobian of the map e tA on K" is equal to e 2rlt ; this is easily deduced from the fact that a divergence 
free flow has Jacobian equal to one ILee031 Proposition 18.18] applied to the map e~ 2t e tA . Also, 
the time t map on R n is one to one and is inverted by e" tA [Lee03| , Theorem 17.8]. It follows that 
the Lebesgue measure of e~ tA B is equal to e~ 2nt times the Lebesgue measure of B. By Lemma 14.11 
the set e~ tA [B] is contained in A n and for each x £ e~ tA [B] the stopping time r x is greater than t 
(apply the lemma for all s < t). The conclusion of the lemma follows. □ 

To see that that yields the special case of Theorem 12.11 let B be a subset of the interior of A„. 
Then the probability, starting with a uniform random point of A n , of being in B at time t is the p n 
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measure of points x € A„ that evolve in time t to a point in B. By Lemma [5TTI this is e~ 2nt p n (B). As 
the probability of a Bin (n, e~ 2t ) random variable being equal to n is precisely e~ 2nt , the conclusion 
of Theorem l2.1l is verified in this case. □ 

The full theorem is a consequence of the next result. Extend the notation for the map 77 so that 
it acts on sets: for B C A m , T C T(t; k ) and S C S(n\ k ) we define 

rj[B, T, S] := {^(y, t, s) : y £ B, t £ T, s £ S} . 

In general, for K C A m x T(t;k) x S(n;k) we denote rj[K] := {^(y, t,s) : (y, t,s) £ K}. 

Lemma 5.2. Let m = n — k < n, lett > 0, and let K = Kq x S(n ; k) be the product of a measurable 
subset of A m x T(t;k) with all of S{n\k). Let g be the density defined in equation (15.21) . Then 

(5.3) Hn{r]{K)) = f g(t) dp L (y) dt. 

J K 

This lemma says that Lebesgue measure on A„ pulled back by g yields the product of normalized 
Lebesgue measure with the measure on T(t\k) having density g , provided that one sums over all 
possible embedding sequences s. Before proving it, let us check that it implies Theorem 12. II Setting 
K = B x T(t; k) makes rj(K) the set of all x for which d> t (x) £ B. The integral in (15.31) factors as a 
product in the y and t variables, resulting in 

hn{K) = Hm(B)p t (n,rn). 

The left-hand side is P(d ) t(x) £ B) and the right-hand side is the normalized Lebesgue measure of 
B times the probability of Bin (n,e~ 2t ) = m. This proves Theorem 12.II 

There are two tricky points in establishing (15.31) . One is seeing how the product of m dimensional 
measure on A m and k dimensional measure on T(t;k ) maps to m + k = n dimensional measure 
on A„. In fact the bifurcation at a variable time s £ [t,t + dt] results in a product with a spatial 
interval proportional to the normal velocity of of the flow toward the boundary of A„. The second 
tricky point is why the total of these factors involving normal velocities is always constant if one 
sums over all sequences s. For the first of these two tricky points we record the following lemma, 
for which all notation is declared to be local. 

Lemma 5.3. Let H be a subset of R d which is locally a halfspace with boundary K. Let B be a Borel 
subset of K, let F be a vector field on H with F ■ n > 0 on K where n is the inward pointing unit 
normal. Let {$ 4 } denote the flow associated with the differential equation x' = Ffx). Let B t C H 
be the set {<I> s (x) : x £ B, 0 < s < t}. Then 
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PROOF: The map $ : K x [0,f] taking (x, s) to < f ) s(x) has Jacobian J(x) = F(x) • n(x) at the point 
(x, 0) because it has a matrix representation in blocks of dimension d — 1 and 1 of 


Id-1 

* 

0 

F(x) _ 


in local coordinates K x [0, oo). The measure of B t is equal to / Bx [ 0t ] J(x, t)(dx x dt). Using 
continuity, J(x, t) = (1 + o(l)) J(x, 0 ) whence the integral is (t + o(t )) f B J(x, 0 ) dx which matches 
the right hand side of (15.41) . □ 


Continuing the proof of Lemma 15.21 we observe that it suffices to show (15.31) for rectangles K = 
BxT. Let B be a small ball about a point y and T be a small rectangle {cj < tj < Cj + A : 1 < j < k} 
about a point t. This makes F ■ n roughly constant over the backward evolution of B jumping at 
times in T, at any fixed time s < t, as long as the embedding sequence a is held fixed. Letting 
tj. |-i = t and to — 0 , we have 


(5.5) 


\v(B,T,a) | 


|£||T|exp 


k 


j)(t j+1 - tj) 

3=0 


k 


i u 


Here, | • | denotes Lebesgue measure, F is the vector field F(x) = A n _j + ix, and A j is F • n evaluated 
at the point where the transition is made from dimension n—j + 1 to dimension n—j. The asymptotic 
equivalence is as the diameter of B goes to zero. The exponential reflects the fact that the backwards 
flow shrinks volume at a rate of — 2 (n — j) over the time interval between tj and t/+i- 


The last step is to sum over er. We do this inductively, beginning with the sum over all values 
of er i. Still supposing that B is a small ball about y and T is a small rectangle about t, we see 
that as time goes backward from t all the way to t \, the backward evolution takes y to a point 
y' = ( y [,... y' n _x) € A„_i. The point y' corresponds to one of n + 1 possible points under the 
different possible choices of o\ = 0 ,..., n, the j th of which is the point (yj,..., y', y',..., y' n _i) 
(where j = 0 or n corresponds to prepending 0 or appending 1). The identification <jj embeds A n _i 
in A„ as the set of points whose j and j + 1 coordinates are equal. The dual vector v i —> v • n is the 
functional v i —> Vj+i — Vj. Therefore, at the point (yj, y ),..., y' n _i), the value A j = —F ■ n is 
(y ' +1 — 2 y'j + y'j) - (y' — 2y' + Vj-i) which simplifies to y ' +1 — y '_ l . Note that we have used —F 
because Lemma |5.3I refers to the velocity of the backward flow. To summarize, 

(5.6) The factor A j is the sum of the two gaps on either side of the coalesced point. 

When j = 0 or n we obtain just one gap. Summing over values of <j\ from 0 to n gives twice the 
sum of all gaps, which is just the constant 2 . 


Inductively, still restricting B to a small ball about y and T to a small rectangle about t, we see 
that the sum over a \ 1 ..., eq of Tl*_i A j is 2*. When i reaches k, we see that the sum over all <r of 
Jd ^ =1 A j is 2 k . We have shown that summing over a leads to a constant factor. Letting E denote 
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the set of all embedding sequences, we see that that ? 7 (y,t,E) := ?7(y> a ) ^ as density 2 k 

independent of y. Integrating over y is therefore trivial and leads, for T a small rectangle, to 


\v{B]T, E)| = |T| • \B\ ■ 2 k exp 


-^2{n- j){t j+1 -tj) 
3 =o 


Recalling that (j,j is j! times Lebesgue measure, we may rewrite this as 

k 


1 


-Mr7(B,T,E)) = \T\ 




ml 


2 k exp 


-tj) 

3=0 


hence 


fJ"n(r](B, T, E)) = 2 k — -\T\fi rn (B) exp 


ml 


~'52 2 (n - j)(t j+ i -tj) 
3=0 


Comparing to (15.211 with k = n — m shows that for small rectangles T, 

Mj7(B,r,E))= fl (t) \T\fim(B). 

Integrating over T yields the result for general sets T, which is (15.31) . implying Lemma 15.21 and 
Theorem 12.11 □ 


6. The time reversed Markov chain on finite configurations 

The evolution of a finite configuration is deterministic. Nevertheless, for some choices of initial 
measure, projecting the configuration to its cardinality produces a continuous time Markov chain. 
The next result is an immediate consequence of Corollarv l2.2l Recall that p n is the uniform measure 
on A n and let A4\ denote a mixture of laws /i n when n has Poisson distribution with mean A. 

Proposition 6.1. Forx £ A, let N(x.,t) denote the number of particles in the configuration at time 
t started from x. 

(i) If x is a random variable with law g, n , then N(x, •) is a pure death Markov chain on Z + with 
initial law 5 n and jump rate at time t from state n given by R(n, t) = 2n (hence transition 
kernel given by (Em;- 

(ii) If x is a Poisson with law M\, then N(x, ■) is a pure death Markov chain on Z + whose initial 
distribution is Poisson(X) and whose transition rate is again 2n. □ 

We now describe the time reversals of these two chains. Both time reversals are pure birth chains. 
The first chain is a time homogeneous chain. It starts at 1 and has rate q(k : t) := 2k of jumping 
from k to k\. It is killed at the time a that it first reaches n + 1. The second has jump rates that 
depend on time but not on state. The initial distribution is Poisson with a specified mean /3 > 0. 
The rate of an upward jump at time s is r(k, s) := 2(3e 2s ds. 
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Proposition 6.2. 

(i) Let v n be the law on right-continuous trajectories of the pure death Markov chain with death 
rate 2k from state k, started from S n and run until the time r that it reaches zero. For a 
trajectory uj, denote the time reversal by s ^ ui'(s) := u>(t — s), made right continuous by 
setting uj'(s) = w'(s + ) when t — s is a jump time of to. Then the v n law of u/ is the same as 
the law of a pure birth Markov chain with birth rate 2k from state k and initial distribution Si, 
on the time interval [0, a), where a is the hitting time on n + 1. 

(ii) Fix t and A > 0 and let Q\ be the law on right-continuous trajectories of the pure death 
Markov chain with death rate 2k from state k, started from a Poisson of mean X and run for 
time precisely t. Let u / again denote the right-continuous reversal of the trajectory u>, that is, 
ui'(s) = ui(t — s) except at jump times. Then the Q\ law of ui' is the same as the law of a pure 
birth Markov chain with initial distribution Poisson of mean (3 := Xe~ 2t and birth rate 2f3e 2Xs 
at time s, run for time precisely t. 

PROOF: Part (i): A trajectory ui is specified by its n jump times, Si, S2, ■ ■ ■, s n - 1, r. The density of 
the trajectory with respect to n-dimensional Lebesgue measure is 

2 "n! exp (— 2 nsi — 2 (n — l)(s 2 — si) — • • • — 2 (r — s n )) . 

Similarly, parametrizing trajectories of the birth chain by the jump times (ri,... r„_i, a) gives a 
density of 

2 ”n! exp (—2 n — 4 r 2 — • • • — 2n(cr — r„_i)) . 

Setting a = r and rj=r — s n -j for 1 < j < n — 1 reduces both exponents to — 2r — 2 Y^j=i s j- 

For part (ii), begin by observing that a death rate of 2 n is equivalent to the n particles each dying 
independently at rate 2. Thus Q\ is the law of Poisson-A many particles each with an independent 
death time whose law has density 2e~ 2s with an atom of size e~ 2t being still alive at time t. A 
trajectory may be specified by death times together with the number remaining alive at time t, 
which is a Poisson process with intensity 2Ae _2s ds + A e~ 2t 6t- The birth chain is specified by these 
same parameters. Clearly the number alive initially is Poisson with the correct mean (3 = Xe~ 2t . The 
arrival process has rate 2/3e 2As = 2Ae -2A ^ -s ) which agrees with the intensity of the point process 
of death times at time t — s. A Poisson process is completely specified by its intensity, establishing 
the distributional identity. □ 

Let R(n , t) be a jump rate for a pure birth process. We define an associated Markov transition 
kernel with state space A as follows. First a sample path of the birth process is generated, giving the 
jump times of the trajectory. Between jumps, the trajectory evolves deterministically via x' = — A„x 
when x E A°. At jump times, a uniform [0,1] random variable U is generated, and the particle 
closest to U splits (including the frozen particles at 0 and 1). Let IP = denote the law of 

a trajectory with jump times given by R, started from p,, with each trajectory a/ stopped at time 
t( W '). 
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Theorem 6.3 (time reversal). 

(i) Let Ri(k,t) := 2k. Define p,* to be the measure giving probability 1/2 to to,o(0) and 1/2 to 
to,i(0). In other words, one of the frozen particles in the empty configuration is chosen by fair 
coin-flip to duplicate. Let r be the hitting time on A n+ \. Then the time reversal of P^t.Ri is 
the law of forward evolutions started from p n and stopped upon hitting Ao. 

(ii) Let i? 2 (/c,t) := 2/3e . Let t be a fixed constant, t. Then the time reversal o/P^^,t,r 2 the 
law of forward evolutions started from A4\ and stopped at time t. 

PROOF: Let IV denote the space of reversals of finite trajectories of forward evolutions. The time 
reversal of the forward evolution from p n stopped at r (respectively the forward evolution from 
stopped at t) is a well defined measure on M\ (respectively M 2 ) on IV. By Proposition 16.21 the 
splitting times for Mi (respectively M 2 ) are given by the birth process R\ (respectively i? 2 )- The 
only thing left to check is the law of the jump at the splitting times. 

By (15.511 . the jumps are chosen independently and proportionally to A j, the normal component of 
F upon embedding by ij. By (15.61) , these are proportional to the sums of consecutive pairs of gaps. 
If 0 = Xq < X\ < • • • Xk < Xk+\ = 1 and U is a uniform random variable on [0,1] then the index of j 
such that Xj is the nearest point to U has law proportional to the sums of pairs of consecutive gaps. 
This completes the proof of the theorem. □ 

7. The infinite system and a weak solution 

In this section we construct the infinite process, which is a probability measure P on the space 
of trajectories of infinite point configurations on R. This requires constructing the space Ai of 
infinite point configurations, then the space of trajectories 12, then giving the probability measure P. 
Although the usual choice for the space of infinite point configurations is that of counting measures, 
we use here a more direct construction that lends itself to taking limits of the finite configurations 
we have already defined. 

Topologize A by treating the projection ir as describing attachment maps from dA n to (Jfc=o ^fc 
for each n. Notice that this identifies all n faces of each dA n . Also note that A is not compact: 
a sequence of points in A° for n —> 00 has no limit. In this topology the trajectories |$t(x)} are 
continuous. They are differentiable except at collision times; one-sided derivatives exist even at 
collision times. 

Let A[, a ’^ denote the simplex of vectors of length n of elements of [a, b] whose coordinates are 
nondecreasing, in other words, A„ with [0,1] replaced by [a, 6]; denote A ^ := A^f L ’ L \ Let A^’ 6 ! 
denote the disjoint union of the interiors, also topologized by attachment. Each A^’ 6 ! embeds 
naturally in A[° ,d l for [ 0 , 6 ] C [c,d\. There are natural projections n [ c ,d],[a, 6 ] : A[ c ’ d l —> A[ a ’ b l that 
ignore points outside of [a, b\ and act as one-sided inverses to the natural embeddings. We denote by 
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7T[a,b] the projection from the inverse limit to Because [a, b\ may now vary, we will sometimes 

need the notation $t,[a,b]( x ) instead of <&i(x) for the time t map on A[ a ’ b l. 

Let Ai denote the inverse limit of these projections as M -A oo. Thus, for us, an infinite point 
configuration 77 is a compatible collection of its finite projections rather than a counting measure. The 
projection 7T[ a ;,] induces a projection from A4 to trajectories on A[ a ’ b l via (7 T[ a ^oj)(t) := 7T[ a ;,](a;(f)). 
Particles in the configuration 7 r[ a . i b]cc(t) die as they pass out of [a, b} and are born as they pass into 
[a,b\; this does not create a discontinuity in because of the identifications t nj o and i n ,n+i 

which allow creation and destruction at the endpoints of [a, b ]. In this context, the Poisson measure 
with intensity A, denoted ), is the probability measure on Ai that projects under each 7 T[ q ;,] to 

v{\ ,M]). 

Let fl denote the space of continuous trajectories on Ai, that is, continuous maps w : M + —> Ai. 
Topologize fl by uniform convergence on compact time intervals [0,T]. Here, by the inverse limit 
construction, uniform convergence means uniform convergence of each projection 7 tlw(-) on each 
compact time interval [0, T], We care only about a tiny subset of fl, namely those trajectories that 
might arise as limits of trajectories $ 7 , 1 , (x). To avoid having to craft arguments for paths that 
might have exotic behaviors, we define this closure, ffo- Formally, for each n and L and each point 
x € A„ , the trajectory {<f> tj L(x)} lifts to an element of Q, the one in 7 rf 1 with no points outside 
of [—L,L\. Let f2° denote the set of such trajectories and let Qq denote the closure of fi° in fi. 

Our first aim is to show that elements of Slo look like trajectories of <I>: they can never uncoalesce 
and they satisfy the evolution rule (ED- There is no probability involed here. A further goal will 
be to see that with probability 1 there are no multiple collisions and that in the limit the number of 
particles remains locally bounded in probability. We begin with formal definitions of collisions, of 
what it means to solve Ob and of how we count the number of particles in a spacetime rectangle. 

Definition 7.1 (collision times). A collision time for a trajectory oj € fl is a time t such that for 
some L,e,n and 1 < j < n — l, the trajectory 7 is in the interior of A^ for the time interval 
(t — e,t) and is in t n — i,jA^f}i for the time interval [t,t + e). The location of the collision at time 
t is the j th coordinate ofoj(t). Although this will turn out to have probability zero , the definition of 
collision allows for more than one collision at time t or for a multiple particle collision. 

Definition 7.2 (solution). 

(i) A trajectory lo £ is said to obey (11.11) if each 7 €E A° obeys (11.11) at time t as long as 
2 < i < n — 1 and t is not a coalescing time for Xi. 

(ii) A weak solution with Poisson-X initial conditions is a probability measure Pa on fl giving 
probability one to the set of trajectories obeying O) and such that the law w( 0 ) under Pa is 
Poisson with mean X. 

(in) A strong solution with Poisson initial conditions is a map S from A i to Q such that 2 ( 7 y)( 0 ) = 77 
and 2 ( 77 ) obeys (11.11) for V(X)-almost every 77 . 
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Definition 7.3 (occupation). Let n(a,t)(u>) = 1 if u>(t) has some coordinate equal to a and zero 
otherwise. Here the trajectory ui may have values in A4 or in any A[ a ’ b l, noting that n(c,t)(n[ ai b] o 
oj) = n(c,t)(ui) when c £ [a, &]. Let n(A,t ) := 'Yh x &A n ^ x '^) denote the occupation of the set A by 
ui(t). Note that for A = {—L, L) the value of n(A , t) is the dimension n of the face A ^ such that 

n L (v(t)) £ A.[ 1 i) . 

Most of the remainder of this section is devoted to proving Theorem 12.31 This will be proved by 
taking a weak limit of solutions to the finite system. 

Proposition 7.4. Paths in Do have the following properties. 

(i) Paths never un-coalesce. Specifically, fix real a < b and 0 < s < t and fix oj £ Do- If 
n({a,b},u) = 0 for all s < u < t, that is, no particle enters or exits [a, b] during the time 
interval [s,i], then n([a, b\, t)(ui) < n([a, 6], s)(w). 

(ii) Paths in Do satisfy (HD- 

Proof: For (i), fix w £ D 0 . Choose a sequence {ujl ■ L £ Z + } converging to oj such that ojl is a tra¬ 
jectory for some point x^ £ A^ L \ Fix s < t and an interval [a, b ] with n({a , b}, u)(uj) = 0 

for all u £ [s, t\. This implies n({a, 6}, u)(ojl) = 0 for L sufficiently large. Denote the initial number 
of particles in [a, b] by r := n([a, b\, s)(w). We need to show that n([a, b\,t ) < r. 

Claim: If x £ A ^ has coordinates Xt = y < z = Xj then the corresponding coordinates of $*,l(x) 
differ by at most (z — y)e 2t . Proof: Up to the first collision time, each gap yk := Xk — Xk -i increases 
at rate at most 2 yk due to the velocity of Xk being at most yk and the velocity of Xk-i being at least 
—yk- This proves the claim up to the first collision time. Induction then proves it for all times. 

Define the £-mesh occupation of [a, b\ of a point x £ A^l, denoted n e (x), to be the least j such 
that the set of coordinates of x can be covered by j intervals of length no more than e. The fact 
that wj —» w implies that for every e > 0 there exists Mo(e) such that 

m>M 0 => n E (wi(s))=r. 

Because the trajectories ujl do not cross the endpoints {a, b} in times in [s, f], each particle at time 
t comes from one or more coalesced particles at time s, hence from the claim we deduce that 

n e 2 (t-. )e .($i ) L(x L )) < n £ (x L ) = r . 

Sending L -A oo and noting that iT[ a ^\CJ is the uniform limit of n a ^\CJL on [s, t], we see that ns(ui(t)) < 
r for any S > e 2 ( t- ^£. Because £ is arbitrary, this implies that n([a, b],t) < r. 

For (ii), let Xi denote a particle for which t is not a coalescing time. Uniform convergence of the 
i—l,i and i +1 components of 4’t,i(xi) implies that both the position and derivative of the particle 
Xi converge uniformly. The time derivative in the limiting trajectory is therefore the limit of the 
derivatives for the finite trajectories, hence the limit trajectory obeys (0 at time t. □ 


COARSENING IN ONE DIMENSION: 


INVARIANT AND ASYMPTOTIC STATES 


15 


Lemma 7.5. 

(i) Let have distribution V(X,L). Any weak limit as L —>■ oo of the laws of the trajectories 
is a weak solution with Poisson-X initial conditions. Consequently if the family of 
pushforwards ofV(X,L) under x i—► {^^(x)} is tight, then such a weak solution exists. 

(ii) If for every t the limit lim^oo ^ t (^L(v)) exists V(X)-almost surely, then this limit defines a 
strong solution with Poisson-X initial conditions. 


PROOF: The weak limit of measures supported on a set S is supported on the closure of S. Therefore, 
any weak limit of the trajectories <& t ,L is in f l 0 and hence, by part (in) of Proposition !?. 4l obeys (11.11) . 
The initial conditions are a weak limit of measures V(X, L) which is the law V(X) on Q. This is all 
that is needed for (i) along with the observation that tightness of this family of laws implies the 
existence of a limit point. Statement (ii) is an immediate consequence of part (Hi) of Proposition !?.!! 
and the definition of a strong solution. □ 


All the time t marginals of the finite system are explicitly known. The only work in establishing 
tightness is to check that the supremum over a time interval s G [ 0 ,t] of the number of particles in a 
fixed interval [a, b] at time s does not have positive mass going to infinity as L —> oo. We have very 
little information about joint distributions of the process at two or more times. However, we can 
accomplish what we need by an identity that bounds the total occupation of the interval in terms 
of particle velocities. 


Proposition 7.6. Let N(a,t)(co) := t] n ( a > s )( w ) denote the cumulative occupation of {a} up 

to time t. Then for any n, any time t and any initial configuration x € 

r L 


-L 


N(x,t)dx < 2 Lt. 


PROOF: Fixing x, the successive collision times ti(x),T 2 (x), ... may be treated as constants. An¬ 
alyzing the trajectories separately on [0, ti), [n, 72 ), etc., and summing the results, we may as¬ 
sume without loss of generality that t < T\. The function I(a,t) is the sum of n indicators of 
graphs of functions, these functions being the trajectories of x ±,..., x n as a function of time. Let 
Vj(t) = 2 Xj(t) — Xj + i(t) — Xj-i(t) denote the velocity of the j th particle at time t, where we have 
used Xj(t) to denote the j th coordinate (<& t (x))j of the system at time t. One has an easy bound 

n n n 

( 7 . 1 ) - X j)+Yl \xj — Xj+i | < 2 L . 

3 = 1 i =1 3 =1 

Let N(a,t,j) = #{s £ [0,f] : Xj(t) = 1} count incidences for just the j th particle. Banach’s 
Indicatrix Theorem | IBan251 Theoremes 1 and 2] (see also INat55| j states that f_ L N(x,t, j) dx is 
equal to the total variation of the trajectory {ar,-(s) : 0 < s < t}. This is just |i>j(s)| ds. Summing 
gives 
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Everything is nonnegative and we may interchange the integral and sum to obtain 

r L rt I n \ rt 


J ^N{x,t)dx = (yjMs)lj ds ^J 0 


2 L ds = 2 Lt. 


Lemma 7.7. For any n, any time t, any L, and any interval [a,b] C [— L, L], 

E f N(x, t) dx < 2(6 — a)t + 2t\~ 1 e 2t , 

J a 

where the expectation is with respect to trajectories whose intial point is uniform on . 


□ 


PROOF: As in the previous proposition, 

N(x,t)dx = f V(s)ds 
Jo 

where V(s) = Ylj |^j(s)| is the sum of the speeds of all particles in [a, b\ at time s. Instead of 
V ( s ) < 2 L for all s, we have 

V(s) < 2(b — a) + (M(s) — b) + (a — m(s)) 

where M(s ) is the location of the first particle to the right of b at time s and m(s) is the location of 
the first particle to the left of a at time s. Each of the terms M(s) — b and a — m(s) is an exponential 
of mean A _1 e 2s , truncated to L — b and a + L respectively. The mean of the truncated exponential 
is at most the mean of the untruncated exponential, leading to 

V{s) < 2(6 — a) + 2A _1 e 2s . 

Plugging this into m proves the lemma. □ 

Lemma 7.8 (velocities are bounded in probability). Fix [a, b] and T. For any trajectory uj := 
{<l) tj i(x)} in n°, let V(ui) denote the maximum absolute velocity of any particle in [a,b\ at any time 
t < T. Then under Poisson initial conditions, V is bounded in probability, meaning that 

V(X,L) ({w : V(w) > M}) < g{M) 

for some function g going to zero as M —> oo and depending on a, b and T but not L. 

PROOF: Label each initial particle Xj by the positive real number £(xj) := e 2T max{x' 7 —Xj-i, Xj +1 — 
Xj}. From the claim proved in Proposition 17.41 we know that at all times between 0 and T, the 
gaps between Xj and Xj -1 is at most £ {xj) provided that we interpret the gap to be zero if the two 
particles have coalesced. The same is true of the gap between Xj and x J+ i . It follows that the two 
gaps adjacent to a particle x at any time up to T is at most the maximum, call it W(x), of £{xj) 
over all particles Xj that have coalesced into x. Because the velocity of any particle at any time is 
bounded by the maximum adjacent gap, the position of any particle x at time t < T must be within 
distance tW (x) of the initial position of some Xj that has coalesced into x. If the particle x at time 
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t < T is in the interval [a, 6 ], and if Xj is the initial position of a particle that has coalesced into x, 
then the distance d(xj, [a, 6 ]) from Xj to the interval [a,b\ must be at most TW{x). Choosing j for 
which W(x) = (.{Xj ), we see that d(xj, [a, 6 ]) < T£(xj). Therefore, 

(7.3) F(w) < 14(w) := T ■ s\xp{i{xj) : d[xj,[a,b\) <Tl{xj)} . 

Couple the distributions of V* under the laws V(X,L) as L varies, by taking the restrictions to 
[—L,L] of a single pick 77 from the Poisson measure 'P(A) on R. The function t{x :j ) for particles of 
ir^r/ are at most what they are for the corresponding particles of 77 , with equality except for the first 
and last particles in [—L,L\. The conclusion of the lemma will therefore follows once we show that 
for infinite configurations, 

V{\) [{77 : V*(V) = 00}] = 0. 

This follows from the fact that the supremum in (17.31) is taken over an almost surely finite set. 
This, in turn, is a simple consequence of the Borel-Cantelli lemma, once one observes that the 
probabilities q n and q' n are summable, where q n is the probability of existence of an initial particle x 
in [b+n,b+n+ 1] with Ti{x) > n and q' n is the analogous probability for a particle in [a — n — 1, a — n]. 
□ 

Proof of Theorem 12.31 Fix L, A > 0. On the probability space of trajectories with initial distri¬ 
bution V{\,L), define random variables F = Y(L,a,b,t) to be the maximum number of particles 
in the interval [a,b\ at any time s E [0, t]. Any particle contributing to this number is either in 
[a — 1, b + 1 ] at time zero or crosses the interval [a — 1, a] or crosses the interval [ 6 , b + 1]. A particle 
crossing an interval J of length 1 contributes at least 1 to f 7 N(x, t ) dx. Therefore, 

P(F > 3 y) < P(F 0 > y) + P(Yi > y) + P(T 2 > y) 

where Fo is the initial number of particles in [a, 6], Yi = J^_ 1 N(x,t) dt and Y 2 = / fe b+1 N(x,t) dt. 
Each of Fo, Fi and F 2 is bounded in expectation by a constant depending on a, b , t and A but not L. 
Therefore, P(F > 3 y) < C/y, finishing the proof of tightness. 

Tightness of these variables for all a , b is tightness of the pushforwards in part (i) of Lemma T7. 51 
the lemma then implying existence of the weak solution which we denote henceforth {di t }. By 
construction the two-sided derivative of each particle’s position exists away from its collision times 
and (ED is satisfied. 

To argue for almost sure distinctness of collision times, fix A and T and the interval [a, b]. For 
e > 0, define the event G e of an e-almost multiple collision by time T on [a, b] to be the event that 
for some times s and t £ [s, s + e] there are distinct pairs of particles whose positions differ by at 
most £ (the two pairs can share one particle but not both). If ujl G f2° with ujl —>• u>, and u> has a 
multiple collision at some time t <T, then for all e > 0 there will be Lq (e) such that L> Lq implies 
ujl has an £-almost multiple collision. Therefore, to finish the proof of the theorem, it suffices to 
show that the probability of G e under any P(A, L) is bounded above by some function n{e) going 
to zero as e —> 0, independent of L. 
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Fix an integer M > 0 and let Hm denote the event that the maximum absolute velocity of any 
particle in [a, 6] up to time T is at most M. Let S be the set of multiples of e in [0, T]. Then, on Hm, 
the event of an e -almost multiple collision implies that for some time s £ S, there are two distinct 
pairs of particles (possibly sharing one particle) within distance (M + l)e at time s. Any pair of 
points in [a, b] within distance (M + l)e of each other is in the same interval [a;, x + 2 (M + l)e] for 
some x = k(M + l)e, 0 < k < 2(6 — a)/{{M + l)e). Therefore, on Hm, the event of an e-almost 
collision is contained in the union of at most |Sj(2(6 — a)e -1 /(M +1)) 2 events that two distinct pairs 
of points at a specific time s are in two distinct intervals [x\, X\ + 2(M + l)e] and [x 2 , £2 + 2(M + l)e], 
together with the |<Sj2(6 — a)e~ l /{M + 1) events that two distinct pairs of points at a specific time 
s are together in the same interval [a?, x + 2 (M + l)er]. 


We bound this from above using the fact that the time s marginal of the number of points in any 
interval of length c is a Poisson of mean A e~ 2s c < Ac. The probability of a Poisson with mean v 
being at least 2 is at most z+/2 and the probability of it being at least 3 is at most u 3 /6. Therefore, 
applying the bound with c = 2(M + l)e, and using the bound |S| < Te” 1 + 1, we see that the 
probability of an e-almost collision is at most 

(Te -i + 1} 2(6 — a) (2(M + l)e) 3 + ^ + 1} (2(b - a)^ 2 f (2(M + 1 )e) 2 ' 2 


M + 1 6 

Clearing away irrelevant stuff, we see that 


2(6 — a) 

£ 


M + 1 






V{\,L)(H M nG e ) < CT(b — a) 2 M 2 £. 

By Lemma 17+1 the probability of H^ is bounded above by g(M). Choosing M = e -1 / 3 we see that 

V{\ L)(G e ) < CT(b - a) 2 e 1/3 + g{e ~ 1/3 ). 


This goes to zero as e 0, finishing the proof of Theorem 12.31 


□ 


8. Further comments and questions 

The most glaring absence of a result concerns strong solutions. The following conjecture, together 
with Lemma 1731 would imply the existence of a strong solution to (ED- A proof seems not too far 
off via results along the lines of Lemma 17.71 

Conjecture 8.1 (strong solution). For P-almost every N £ fl, the limit 

V t (N) := hm V { t L) (N) 

L—y oo 

exists and defines a trajectory t i-A H >t{N). 

Conjecture 18.11 also implies the following restatement of Theorem 12.41 in terms of the evolution 

W- 


(8.4) 


V(X) o T" 1 = V(Xe~ 2t ). 
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Numerical simulation suggests a substantially more general result than those obtained in this 
paper. We believe that the Poisson distribution is an attracting fixed point for the dynamics: any 
reasonable initial measure, if allowed to evolve under these dynamics and rescaled, should converge 
Poisson. There is some numerical evidence for this. In one experiment, 250 million points were 
initially placed on the unit interval with uniform density and periodic boundary conditions. New 
points were then placed halfway between pairs of adjacent points, and the initial points were removed. 
The resulting configuration is known as the ID Poisson-Voronoi configuration. Its cell sizes are one- 
dependent but not independent, with marginal density 4x exp(— 2x). Figure [8] illustrates histograms 
of the normalized cell sizes as time grows. As increasingly many of the cells disappear, the pdf 
approaches the exponential distribution, the cell size marginal for the Poisson. The data is also 
compatible with asymptotic independence of neighboring cell sizes. Additional supporting data 
from this initial configuration and several other initial configurations are reported in |Lazll| . 



To make this into a precise conjecture, let v be a probability measure on R with finite mean and 
let Q u be the stationary renewal process on R with renewal distribution v. Intuitively, Q v is the law 
of a random configuration of points whose gaps are IID u. 

Conjecture 8.2. For any v on R with mean m < oo, there is a weak solution to with 

initial conditions Q v . Let vt be the time t law of such a weak solution, with space rescaled by e~ 2t . 
Then Vt —>• V(l/m) weakly as t —>• oo. 


This conjecture is the concrete assertion of the universality mentioned in the abstract. While 
a proof seems farther off for this conjecture than for Conjecture 18.11 we hope to attack it via the 
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dual process introduced in Proposition 16.21 While the forward evolution is deterministic, the time- 
reversed process has randomness. We believe this will help us to show that running the dual process 
back in time t units produces, as t —> oo, a configuration asymptotically independent from any 
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